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Abstract 

Wu and Sprung {Phys. Rev. E, 48, 2595 (1993)) reproduced the first 500 nontrivial Riemann 
zeros, using a one-dimensional local potential model. They concluded — and similarly van Zyl 
and Hutchinson {Phys. Rev. E, 67, 066211 (2003)) — that the potential possesses a fractal 
structure of dimension d = |. We model the nonsmooth fluctuating part of the potential by 
the alternating-sign sine series fractal of Berry and Lewis A{x,j). Setting d = |, we estimate 
the frequency parameter (7), plus an overall scaling parameter (o") we introduce. We search 
for that pair of parameters (7,c) which minimizes the least-squares fit S'„(7, o") of the lowest 
n eigenvalues — obtained by solving the one-dimensional stationary (non-fractal) Schrodinger 
equation with the trial potential (smooth plus nonsmooth parts) — to the lowest n Riemann 
zeros for n = 25. For the additional cases we study, n = 50 and 75, we simply set a = 1. The 
fits obtained are compared to those gotten by using just the smooth part of the Wu-Sprung 
potential without any fractal supplementation. Some limited improvement — 5.7261 vs. 6.39207 
(n = 25), 11.2672 vs. 11.7002 (n = 50) and 16.3119 vs. 16.6809 (n = 75) — is found in our 
(non-optimized, computationally-bound) search procedures. The improvements are relatively 
strong in the vicinities of 7 = 3 and (its square) 9. Further, we extend the Wu-Sprung semiclassical 
framework to include higher-order corrections from the Riemann-von Mangoldt formula (beyond 
the leading, dominant term) into the smooth potential. 
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I. INTRODUCTION 



In summarizing the results of their paper, "Riemann zeros and a fractal potential" , Wu 
and Sprung stated that "we have found analytically a one-dimensional local potential which 
generates the smooth average level density obeyed by the Riemann zeros. We have then 
shown how any finite number of low lying Riemann zeros can be reproduced by introducing 
fluctuations on top of the potential. The mystery of how a one-dimensional integrable system 
can produce a 'chaotic' spectrum is resolved by adopting the concept of a fractal potential 
which, in the infinite limit, would lead to the system having a dimension larger than 
one" [H p. 2597] (cf. |3l IH El El El E])- ("Indeed . . . finding an Hermitian operator whose 
eigenvalues are [the Riemann zeros] may be impossible without introducing chaotic systems" 
p. 3].) 

The Wu-Sprung potential V — which generates the smooth average level density obeyed 
by the Riemann zeros — satisfied Abel's integral equation p^J eq. (6)], and was written 
implicitly as [I] eq. (7)] (cf. [H] [IHl sec. 4]), 

Here Vq = 3.100737r ^ 9.74123. 

Our objective is to reproduce, as best we can, the fluctuations on top of the potential 
Vws{^)j implicitly given by Q, so that the application of the Schrodinger equation to the 
so-amended (smooth plus fractal) potential would yield the Riemann zeros themselves. For 
our exploratory purposes, we adopt (being a particular case of a deterministic Weierstrass- 
Mandelbrot [WM] fractal function) the alternating-sign sine series of Berry and Lewis [TTl 
eq. (5)], 

(—l)"^ sin 

A{x,^) = E^^_ J , (l<d<2,l<7). (2) 

Here, d is the fractal dimension, which — following the box-counting argument of Wu and 
Sprung [1] (cf. [I2]) — we take to be |. We have, in this = | Berry-Lewis context, a 
specific case, 

^(7x,7) = -7^A(a;,7), (3) 

of the "affine scaling law" pTI eq. (3)]. We also scale — in the first {n = 25) of our three 
sets of analyses (n = 25, 50, 75) — A{x, 7) by a parameter a, where n is the number of the 



lowest Riemann zeros we aspire to fit (sec. II A). For the cases n = 50 (sec. II B) and 75 
(sec. II C), we will simply set cr = 1. In sec. Ill, we demonstrate how to incorporate more 

for the cumulated number of Riemann 



terms of the Riemann-von Mangoldt formula 
zeros than Wu-Sprung themselves did, using a semiclassical argument, in deriving xwsiV)- 
(It remains, however, to numerically implement these last findings.) 



II. ANALYSES 



A. n = 25 



We proceed, to begin, trying to fit the first twenty-five Riemann zeros by finding dis- 
tinguished values of the two parameters (7 and a). We randomly generate trial values 
1 < 7 < 10 and < a < 10. (Numerically-speaking, we truncate the summation in ^ by 
summing from m = —30 to m = 30 (cf. [HI App.]). We have not yet gauged the sensitivity 
of the various results in this paper to this choice of cutoff — nor to the setting = | nor 
further to the specific measure of fit (sum-of-square deviations) employed — though it would 
certainly be of interest to do so for any or all of them.) 

If we use the smooth potential given by ([T]) itself — without any fractal supplementation 
— the sum-of-squares deviation of the first twenty-five eigenvalues yielded by application 
of the one-dimensional stationary Schrodinger equation from the first twenty-five Riemann 
zeros is 6.23907 (Fig. [T]). (This is only 0.0069 percent of the total [non-fitted] sum of squares 
of the zeros themselves, that is, 92569.63, so one might aver that the semiclassically-based 
smooth Wu-Sprung potential is notably successful in well-approximating the Riemann zeros. 
It is, of course, our objective here to reduce this small percentage even further. Let us also 
note that a referee suggested that the scatter in Figs. [T][7] and |9] might be reduced if the 
modulus of the scatter were to be plotted.) 

We randomly generated 4,007 pairs of (7, a) from the indicated ranges, and solved for 
each such pair, the Schrodinger equation with the smooth potential plus the fractal aA{x, 7) 
with the corresponding choice of 7 and a. Further, we calculated the corresponding sum- 
of-squares (5*25(7,0")) deviations from the first twenty-five Riemann zeros. We obtained a 
range of values 5*25(7, ^ [5.7261, 166.075]. In Fig. [2| we display the results for those 1,833 
pairs of the 4,007 that yielded 5*25(7,0") < 15. 
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FIG. 1: Deviations from the first twenty- five Riemann zeros of the first twenty- five eigenvalues 
obtained by solving the Schrodinger equation using the (smooth/non- fractal) Wu-Sprung potential 
Q- The sum-of-squares of these deviations is 6.39207, while the sum-of-squares of these Riemann 
zeros is much larger, 92569.63. 




FIG. 2: Scatterplot showing the sum-of-squares goodness-of-fit statistic 5*25(7,0") as a function 
of the frequency parameter (7) and the scaling parameter (cr). Those 1,833 points of the 4,007 
sampled for which S'25(7,o") < 15 are included. 

Additionally, in Fig.[3]we show the results for those 210 pairs yielding 5*25(7, < 6.39207 
— that is, those pairs which yield results superior (numerically inferior, that is) to those 
obtained with the original smooth unsupplemented potential ([T]). Fig. |4]is a histogram of the 
values of 7 — having been uniformly sampled from [1,10] — occurring in these 210 pairs. 
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FIG. 3: Truncation of the previous scatterplot to those 210 (fit-improving) pairs yielding 
5*25(7, 0-) < 6.39207. For a// such pairs a < 2.86475. 



# 




FIG. 4: Histogram for the frequency parameter 7 corresponding to those 210 pairs (7, a) for which 
5*25(7,0') < 6.39207, the value obtained from solving the Schrodinger equation using the smooth 
unsupplemented Wu-Sprung potential. The classification bins for 7 G [1, 9.5] are of length 

(There were no 7's recorded for a number of bins, including [9.5,10].) The two topmost peaks 
(corresponding to the classification bins 7 G [2.5, 3] and [8.5,9]) may possibly reflect (cf. ([3])) 
an assertion of Berry and Lewis, made (using t for what we denote by x) in regard to what 
they termed deterministic Weierstrass-Mandelbrot fractal functions {W{t)). They note that 
"the whole function W can be reconstructed from its values in the range to < t < 7^0; for 
example, W in the ranges 7^0 < ^ < 7^^o aiid 7~^to t < are magnified and diminished 
versions, respectively, of W in the range < t < 7^0 • • • The repetition and resolution of 
features at to,7^o,7^^o etc, is again obvious" [HI p. 461]. (Similar 7-histograms — Figs. 
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and 17 — will be obtained for the n = 50 and 75 analyses further below. 



In Fig. [5] we have an analogous histogram based on the scaling parameter a. (For no 
value a > 2.86475, though we sample uniformly from [0,10], do we obtain an improvement 
by using aA{x, 7).) In Fig. |6] we have — for these same 210 pairs — a plot of 7 vs. a. (The 
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FIG. 5: Histogram for the scaling parameter a corresponding to those 210 pairs (7,c) for which 
S2^{l,a) < 6.39207 
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FIG. 6: Plot of those pairs (7,0") for which 5'25(7,cj) < 6.39207. The correlation coefficient is 
-0.264263. 

associated correlation coefficient is negative, that is, -0.275492.) 

The minimum over all the 4,007 pairs was S'25(l. 54523, 1.95798) = 5.7261. (The next 
smaller value was S'25(l. 15274, 1.57931) = 5.81754. All other sum-of-squares deviations 
exceeded 6.0.) In Fig. [T] (cf. Fig. [T]) we show the deviations of the predicted eigenvalues at 
this point (1.15274,1.57931) from the corresponding Riemann zeros. In Fig. |8] we display, 
for this same minimizing point, the corresponding twenty-five eigenfunctions drawn at the 
twenty-five eigenvalues for the associated (smooth plus fractal) potential. 



B. 



n = 50 



Now, we present in Fig. |9] the extension of Fig. [T] from the first twenty-five to the first 
fifty Riemann zeros. (It is considerably more demanding to solve the Schrodinger equation 
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FIG. 7: Deviations from the first twenty-five Riemann zeros of the first twenty-five eigenvalues 
obtained by solving the Schrodinger equation using the (non-fractal) Wu-Sprung potential ([T]) 
supplemented by the scaled Berry-Lewis fractal function 1.95798j4(x, 1.54523). The sum-of-squares 
of these deviations is 5.7261, reduced from 6.39207 for the non-supplemented smooth potential 
(Fig.[l]). 
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FIG. 8: Eigenfunctions drawn at the eigenvalues obtained by solving the Schrodinger equation for 
the potential which yields the minimum 525(1.54523, 1.95798) = 5.7261 of the 4,007 pairs sampled 



for this increased n.) The associated sum-of-squares is 11.7002, which is but 0.00261 percent 
of the total sum-of-squares, 448704.56, of the first fifty Riemann zeros themselves. 



In Fig. 10 we show the results of a "rankit" test (described in the Wikipedia on-line 
encyclopaedia) for the normality of the distribution of deviations in Fig. [9] If the observations 
do come from a normal/Gaussian distribution, we expect the plot (necessarily non-decreasing 
in any case, Gaussian or otherwise) to be a straight line. So, there appears to be some 
deviations from strict normality, particularly in the tails of the distribution (suggesting that 
perhaps sums-of- squares might not be the most robust measure of deviation to be employed 
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FIG. 9: Deviations from the first fifty Riemann zeros of the first fifty eigenvalues obtained by solving 
the Schrodinger equation using the (non-fractal) Wu-Sprung potential ([T|. The sum-of-squares of 
these deviations is 11.7002. 
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FIG. 10: A "rankit" plot to assess the possible normality of the distribution of deviations of the 
fifty eigenvalues from the Schrodinger equation solution using the smooth Wu-Sprung potential 
from the Riemann zeros themselves. For normally-distributed observations, one expects a straight 
line. 

for evaluation of fits — though it is certainly the most conventional and familiar measure). 



We repeated for the case n = 50, the form of analysis conducted for n = 25 (sec. II A), 
but now omitting (in the interest of interpretational simplicity) the scaling parameter (cr) 
— effectively setting it to unity. Based on 1,013 random choices of 7 lying in [1,10], we 



obtained Fig. 11 The minimum achieved was 5*50(3.10007, 1) = 11.2672, while without any 



fractal supplementation at all, the (larger) value of 11.7002 was obtained. (As a matter of 
curiosity, we found that ^50(1, ^a«) = 11.7641 > 11.7002, where qau = ~ 1.61803 is the 



golden mean [14J. Interpolation by third-degree polynomials of the points in Fig. 11 suggests 
that the actual minimum of 11.2556 would be achieved for 7 = 3.092.) The penultimate 
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FIG. 11: Scatterplot showing the sum-of-squares goodness-of-fit statistic 550(7, l)i based on 1,013 
randomly chosen points < 7 < 10. Points below the 7-axis correspond to improvements in fit. 

minimum, 5*50(3.10051,1) = 11.2685, was nearby. Again (cf. Fig. |4]), we appear to discern 
a manifestation of relative minima in the vicinity of both 7 = 3 and (its square) 9. (Also 
around 7 = 6, as in Fig. |4] too. The associated relative minimum of 11.678 occurs at 
7 = 5.96158.) The relative minimum (11.4906) in the vicinity of 7 = 9 was attained at 
7 = 8.78602, so, at that point, one has ^/7 = 2.96412. For the smaller values of 7, the 
plot is very scattered, but for the higher ones, a well-defined curve emerges. (The overall 
maximum is 12.1412 at 7 = 1.8165, while the maximum in the upper range [8, 10] is 12.108 



at 7 = 9.7516.) We apphed (third-order) interpolation to the data in Fig. 11 and obtained 
for the most salient domains of improved fit, 7 G [3.03205,3.1537] (containing our overall 
minimum) and [3.33583,3.38488], [5.85954,6.08565] and [8.45962,9.12272] (containing our 
three distinct relative minima [Fig. ITT]). 



In Fig. 12 (cf. Fig. H| we show the histogram based on those (fit-improving) 219 of 



the 1,013 values of 7 for which S'5o(7, 1) < 11.7002. In Fig. 13 we plot the deviations 
(Asmooth) from the first fifty Riemann zeros obtained using the smooth potential against those 
deviations {A fractal) using the supplemented potential which minimizes the fit, obtaining 
basically a linear relationship. In Fig. [14] we have the counterpart of Fig. [8| for the case 
n = 50. 



C. n = 75 



Fig. 15 is the further extension of Figs. [T] and [9] to the n = 75 case. The associated 
sum-of-squares is 16.6809, which is but 0.0014308 percent of the total sum-of-squares of the 
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FIG. 12: Histogram corresponding to those 219 values of the frequency parameter 7 for which 
5*50(7, 1) < 11.7002, the value obtained from solving the Schrodinger equation using the smooth 
unsupplemented Wu-Sprung potential. The classification bins for 7 G [1, 9.5] are of length |. 
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FIG. 13: The predicted eigenvalues minus the first fifty Riemann zeros: the horizontal axis based 
on the smooth potential and the vertical axis based on the sum-of-squares minimizing fit 




FIG. 14: Eigenfunctions drawn at the eigenvalues obtained by solving the Schrodinger equation 
for the potential which yields the minimum 5*50(3.10007, 1) = 11.2672 of the 1,013 points sampled 
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FIG. 15: Deviations from the first seventy-five Riemann zeros of the first seventy-five eigenvalues 
obtained by solving the Schrodinger equation using the (non-fractal) Wu-Sprung potential ([T]). The 
sum-of-squares of these deviations is 16.6809. 
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FIG. 16: Scatterplot showing the sum-of-squares goodness-of-fit statistic 575(7,1), based on 853 
randomly chosen points < 7 < 10. Points — of which there are 172 — below the 7-axis correspond 
to improvements in fit. 

first seventy-five Riemann zeros, 1.1658469 ■ 10^. We generated 853 values of the frequency 
parameter 7, uniformly sampling from the interval [1,10]. The best fit of 16.3119 was 
achieved for 7 = 3.1106. (The penultimate minimum of 16.3427 was nearby at 7 = 3.11632, 



while the overall maximum [worse fit] was 19.27 at 7 = 2.47689.) In Fig. 16 we plot these 
853 values as a function of 7. (In the vicinity of 7 = 9, the relative minimum of 16.5983 is 



achieved at 7 = 8.64589 ~ 2.94039^.) In Fig. 17 we show the corresponding histogram. 

We continue to add randomly generated points to this n = 75 analysis, and may possibly 
undertake an n = 100 study too. 
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FIG. 17: Histogram corresponding to those 172 values of the frequency parameter 7 for which 
5*75(7, 1) < 16.6809, the value obtained from solving the Schrodinger equation using the smooth 
unsupplemented Wu-Sprung potential. The classification bins for 7 G [1, 10] are of length 1. 

III. HIGHER-ORDER CORRECTIONS TO THE WU-SPRUNG POTENTIAL 



In their semiclassical analysis, Wu and Sprung [T] took into account only the leading term 
of the Riemann-von Mangoldt formula [13], but nevertheless, as our results indicate, doing so 
is able to yield, via the Schrodinger equation, eigenvalues that quite closely approximate the 
Riemann zeros themselves. It might prove beneficial, in trying to account for the (relatively 
small) residual variation — especially since we have been working in a non-asymptotic regime 
— to extend the Wu-Sprung approach by incorporating the higher-order non-oscillatory 
terms of the Riemann-von Mangoldt formula too (cf. [TH] eq. (5)]). 

The (smooth) leading term — of the Riemann-von Mangoldt formula for the number of 
zeros below E — employed by Wu and Sprung took the form 

7 



E E 
N(E) = — log 



(4) 



(This formula pertains to the important [Connes/Berry-Keating — absorption/emission spec- 
trum] "sign" problem [TTl eqs. (6), (11)].) The remaining part (having both smooth and 
nonsmooth components) not utilized by Wu and Sprung is expressible as [131 P- 436] 



where 



S{E) + -5it), 
n 



S{E) = ^^Tgal + ^E), 



(5) 



(6) 
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is the argument function and ({s) is the Riemann zeta function. Further, 

and p{u) = 1/2 — {m}. (Here {u} is the fractional part of u.) The function S{E) is itself 
strongly oscillatory (changing sign an infinite number of times). 

Incorporating the two additional (higher-order — 0{E~'^) and 0{E^^), respectively) non- 
oscillatory (monotonically decreasing) terms of the Riemann- von Mangoldt formula [13] , one 
could seek to find the potential based on 

iV(E) + |log(l + ^)+iarctan^. (8) 

The corresponding potential can, in fact, be straightforwardly constructed (in implicit form), 
using the formula for Abel's integral equation [TH], following the procedure outlined in [T]. 
It takes the form 

xwshK^) 2 2 ^ 2 y^) 

^tanh-1 ( l^tanh-i ( 



I 



log (l + ^) + vr (log (47r2) - 2 log ^V^o 



271 



2 



„• ]„„ ( {l~i)Vo+y/V , 2iv/F^ 

> ((l+i)v^+j)V2j-4Vb 4VV+(2+2i) 



2 log 

+ 



27r V2« - 4Fo 

(-i+t)Vb-tv^ 2vy-yo 



((l+i)v^-l)V-4Vo-2i 4v^-(2-2i) 



27rV-4K) - 2i 

, (l-i)(Vi=^v^F^-V2yo)+V2v^ 
^ > (2i-(2+2i)v^)Vi=2W 



+ 



27r V2« - 4\/o 

•1 ( »((l+OVb+\^) (1-»)VF^ 

^ °S > ((l+j)v/y+l)V-4K,-2i (2+2j)VV+2 

27rV-4K) 
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We have not so far been able to numerically invert xvysnl^), but in Figs. 18 and 19 



we show the difference Xwsiy) — ^wsni^) ^^e ratio 



original Wu-Sprung analysis [T], we set Vq = 9.74123 in both cases.) 



, respectively. (As in the 
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FIG. 18: The original (implicitly expressed) Wu-Sprung potential ([T| minus the one Q with 
higher-order Riemann-von Mangoldt corrections 
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FIG. 19: The original (implicitly expressed) Wu-Sprung potential ([T]) divided by the one Q with 
higher-order Riemann-von Mangoldt corrections 



IV. CONCLUDING REMARKS 

It might be possibly worthwhile exploring, in our context, the use of fractals other than 
the specific Berry-Lewis alternating-sign sine series g (cf. Ha Chap. 12] 1201 ED). (It 
appears that A{x,'y) increases with x, while the Wu-Sprung plots of the fluctuations of the 
fractal potential from the smooth one [1] Fig. 2] do not seem to exhibit such a growth 
phenomenon (cf. [191 Fig- 12.13]). However, certain quite preliminary efforts of ours 
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to explore along such lines have encountered some so-far not well-understood numerical 
difficulties.) Along with the alternating-sign sine series {A{x,'~f)) we have employed, Berry 
and Lewis too analyzed a certain companion (also deterministic Weierstrass-Mandelbrot) 
cosine series C{x,'j). However, this function can only assume nonnegative values, so — 
unless translated — it does not provide an immediately suitable model of the Wu-Sprung 
{both negative and positive) fluctuations. 

In the asymptotic limit, the Wu-Sprung potential Vwsi^) behaves as P, eq. 

vr^x^ r / /tt |x I \ 1 
Vws..y^,{x) = — [LW(^^-— y , 

where LW(-) represents the Lambert-W function [22j. It would be interesting to explore 



|x| — > oo, (10) 



the question of whether the Schrodinger equation can be exactly solved with (10) as a 
potential (cf. [23]). (In regard to this matter, M. Trott commented that "I would guess 
it is enough to look for the asymptotic region. If an analytic solution exists there, one 
might be able to 'guess' the correct one in terms of Lambert-W. If one can't find an exact 
solution asymptotically, then probably no exact one in Lambert-W does exist either".) If 
so, then with the use of perturbation theory [21] , one might possibly be able to expedite the 
computational procedure we have employed above (in which we have solved the Schrodinger 
equation — a demanding task — ab initio for each new set of random parameters). 

Questions pertaining to the nature of the eigenfunctions obtained and of the importance 
of tunneling contributions in the fractal potential remain to be addressed (cf. Figs. M and 



14|). 

Castro and Mahecha — in their proposed supersymmetric model of the Riemann ze- 
ros — suggested the use of a ''fractal [emphasis added] SUSY QM equation instead of 
factoring the ordinary [emphasis added] Schrodinger equation studied by Wu-Sprung" 
[20] p. 788] (cf. [211 ESI EHj). They further proposed the use of a fully general 
Weierstrass-Mandelbrot fractal to model fluctuations in the supersymmetric potential. The 
alternating-sign sine series fractal ^ we have employed in this study is a specific in- 
stance of the WM-fractal [TT], involving far fewer parameters — which is largely why 
we have employed it here in our exploratory numerical analyses. (References to further 
related work of Castro can be found at the number theory and physics archive website 



http: / / secamlocal.ex.ac.uk/people/staff/mrwatkin/zeta/physics.htm[) 



The analytical approach of Wu and Sprung [T] , which has been the basis of our study here 
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incorporates the familiar ("Berry-Keating" |27]) form of density-of-states N{E) given in (|4]). 
However, in Connes's adelic (absorption spectrum) approach the density-of-states takes 
another (cutoff (A)-dependent) form [T71 eq. (11)], 

NconnesiE) = £ log _ ^ log ^ - l) . (11) 

The question, then, arises of whether or not the Wu-Sprung framework can be meaning- 
fully/differently recast in the Connes adehc setting. (Berry has been quoted to the effect 
that "we think this is a matter of formahsm and everything you can write as an absorption 
you can write as an emission spectrum if you manipulate things the right way" [291 p. 247]. It 
has been proposed that the approach of Connes can be implemented using pseudo-Hermitian 
Hamiltonians for which the eigenvectors are null [30]. Shudo et al [31] measured absorption 
spectra and investigated the distribution of eigenfrequencies in "L-shaped" resonators.) 

In the context of fractal strings and sprays, the classical Riemann zeta function "can be 
viewed as the spectral zeta function of the unit interval" [32, p. 2]. 

In their study of the Weierstrass-Mandelbrot fractal function, upon which we have 
strongly relied. Berry and Lewis found, after some delicate analysis, that an attractive 
inverse-square potential generates the "Weierstrass spectrum" [TH sec. 5]. They further 
speculated that there might be no other such potential. 
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